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INTRODUCTION. 

The object of the experiments described in this report, submitted to the National Advisory 
Committee for Aeronautics for publication, was to compare the damping coefficients of an air- 
foil as calculated from a knowledge of the static characteristics of the section with those obtained 
experimentally with an oscillator. The damping coefficients so obtained, according to the con- 
ventional notation, can be considered either as due to pitching or due to yawing, the oscil- 
lator in these experiments being so arranged that the surface oscillates about a vertical axis. 
This is in reality the case when the airplane is yawing about the standard Z-axis, but it can also 
be considered as a pitching motion when the model is so rigged that its standard Y-axis becomes 
vertical. This horizontal oscillation has the advantage of eliminating the gravity action and 
avoiding the use of counterweights, whose presence in the wind tunnel is undesirable because 
of their interference with the air flow. The experimental work was all done in the four-foot 
wind tunnel at the Massachusetts Institute of Technology, in connection with the preparation 
by the writer of a thesis in the course in aeronautical engineering at that institution. 

The apparatus used in the experiments is essentially a bifilar suspension. In designing this 
apparatus great difficulty was encountered in making it so as to obtain moment of inertia 
large enough to keep the system oscillating for a time long enough to be measured accurately 
and without elaborate special apparatus. After several trials it was found necessary to employ 
adjustable counterweights oujpde the tunnel. 

The real point of these experiments was to separate the damping due to rotation from that due 
to translation. Consider the motion of a surface situatedat a distance behind the center of rotation 
about which it oscillates. That motion can be considered as the resultant of two component 
motions: First, a rotation about the center of pressure of the surface; second, a translation 
perpendicular to the wind direction. The larger part of the damping is due to the translational 
motion, and only this part can readily be calculated from static tests. By varying the distance 
between the center of pressure and the center of rotation on the oscillator, the variation of 
damping moment can be observed, and the rotational and translational effects can be sepa- 
rately determined. 

SUMMARY. 

In the first part of the present work, dealing with theoretical damping coefficients, a brief 
discussion of the method of calculation is given. Owing to the limited amount of time no attempt 
was made to test a large series of models. An 8 by 2 inch flat plate was first tested, followed by 
a tail piece with two elevator settings. Static tests were first performed on these surfaces for 
three speeds, 30, 20, and 10 miles per hour, and the corresponding damping coefficients calcu- 
lated. Owing to the inaccuracy of force measurements and the slow damping at very low speeds 
the 10-mile runs for the tailpiece were omitted. 

Dynamic tests were made and experimental- damping coefficients were calculated for the 
same settings used in the statical tests. These tests are included in Part II of the thesis. A 
comparison of results from Parts I and II is given in Part III with brief discussion and conclusion. 

The results were such as to encourage the continued use of the conventional method of 
calculation of damping coefficients. The agreement between the experimental and calculated 
values was extremely good for all positions which would be likely to be occupied by an actual 
tail surface, the maximum difference for such positions being less than 10 per cent. In general, 
although not in all cases, the experimental value was a very little below the calculated. 
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PART L 

THEORETICAL DAMPING COEFFICIENTS. 

Consider a tail surface situated at a distance behind the center of gravity of the airplane 
and rotating about the C. G. If the machine has a forward velocity V and an angular pitching 
velocity 2, the resultant wind velocity relative to the tail surface is V r as shown in the figure. 

The angle between V T and the direction of flight, measured in radians, is tan or tan -p X57.3 

in degrees. 

The change in lift due to the rotation is: 

-SO")*** 

where a is the angle of attack of the tail surface and S is its area. 

Moment- -^|x67.3x^X jS V*Xl 
--^X57.3 fSxPxVxq 

Therefore ^ ^X57.3±SPV 

" ~ 57 - 3 ^ x T ft " fcs./radian/sec. 

where ^ is the slope of the lift curve plotted against angle of attack in degrees. The value of 

^ is very nearly constant in the neighborhood of zero angle, as the lift curve there is approxi- 
mately a straight line. With this value obtained from the statical tests of the section, the 
theoretical damping coefficients due to the translational component can be calculated. 

The statical tests on the flat plate were made at three speeds: 30, 20, and 10 miles per hour. 
A complete set of characteristic curves are given for 30 miles per hour, but only the lift curves 

are given for the 10 and 20 mile runs, as required for the purpose of computing The value ^ 

of course varies very closely as the square of the speed. 

The mean position of center of pressure for the fiat plate is taken at one-fourth of the chord 
from the leading edge — that is, 0.5 inch. The mean position of center of pressure for the tail- 
piece with — 30° elevator setting is taken at the hinge, and that for the zero degree setting is 
taken 0.3 inch forward of the hinge. With these mean positions of center of pressure the value 
of 7. is properly taken and the theoretical value of damping coefficients can be calculated. The 
results are tabulated in Table I and plotted in figures 1, 2, and 3. 
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STATICAL TESTS ON PLAT PLATE AND TAILPIECE. 

The tests were made in the usual manner on the wind tunnel balance, lift, drag, and moment 
readings being taken. The results for the tail surfaces are fully shown in the curves which 
follow. 

I. 

Theoretical damping coefficients of tailpiece. 



[Elevator setting -30 and Q*. Wind speeds 20 and 80 miles per hour.] 



Elevator angle. 


-30'. 


-80". 


o\ 


0". 
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M. 

-v- •••••• 
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part n. 

EXPERIMENTAL DAMPING COEFFICIENTS. 

METHOD OF EXPERIMENTING AND DESCRIPTION OF APPARATUS. 

The general scheme employed in testing was to mount the model in such a "way so that it 
oscillated about a fixed vertical axis. As mentioned in the introduction, the use of a vertical 
axis has the advantage of avoiding the action of gravity and the use of counterweights. For 
minimum interference, a minimum of frictions! damping, and simplicity of contraction a 
bifilar suspension was chosen. 

Referring to the accompanying drawing (fig. 7), it is seen that two beams shaped to 
stream-line form constitute the principal structural members of the set-up. These beams are 
secured to the channel wall by inserting the tongues on the ends of the beams into the sheet- 
metal sockets on the wall, thus making it possible to erect or dismount the apparatus quickly. 
Two piano wires with one end hooked to the upper beam (1), passing through the hole bored in 
the lower beam and running oven the pulleys (9) , are led to the outside of the channel. It is the 
chief aim of the design of the apparatus to put as few as possible of the parts inside the tunnel. 
Two weights of about 10 pounds each are attached to the other ends of the wires. These two 
stretched wires therefore function as a spring under constant tension which supplies the restoring 
moment necessary to keep the model in regular periodic motion. 

The turnbuckles (11) can be adjusted so that the wires pass through the small holes on 
the floor of the channel freely without touching the sides of the holes. 

The oscillating bar (3) is made of two semicircular rods clamped together by screws. The 
bar is clamped to the piano wires by screws at one end, and the model is clamped between the 
halves at the other end. The method of mounting the model is illustrated in the drawing. 

The damping and statical moments in these experiments were both unusually large on account 
of the long lever arm (in the extreme case more than 15 inches). This necessitated the use of 
heavy inertia weights. As shown in the drawing, these weights are placed outside the tunnel 
and are clamped to the rod (18), which is connected rigidly to the oscillating bar (3) through 
the vertical rod (13). At the intersection of the vertical rod 13 and the oscillating bar 3 a 
universal joint is provided, so that the stretched wires and the rod impose no lateral constraint 
on each other and there is no side thrust on the bearing at the bottom of the rod except that 
due to friction in the universal and to the air resistance of the rod itself. 

At the lower end of the vertical rod 13 a pivot rests in a socket. This bearing is intended 
only for constraining the motion of the vertical rod and not for taking any load, as virtually 
all the weight is carried by the piano wires. 

A pointer is attached to one end of the rod (18), and the angular displacements are read 
directly on the dial (20), which is graduated in degrees. 

The apparatus was set up about 4 feet upwind from the balance, so that the stand (23) was 
placed very near to the motor rheostat. Only one observer was thus needed for the experiment, 
regulating the wind speed and taking the reading at the same time. 

SOLUTION OF THE DIFFERENTIAL EQUATION OF MOTION. 

The mathematical principles involved in these experiments are simple and need only be 
summarized. 

The damped harmonic motion. — In this oscillating system the damping is due to two causes, 
the mechanical friction of the mechanism and the action of the air on the system. The resisting 
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moments due to these two causes, together with restoring moment of the piano wires at any- 
instant of the motion, must be equal to the product of the mass moment of inertia by the 
angular acceleration: 

J~at> — m ~ l Tt a) 

where — is the moment of inertia of the entire oscillating system measured in slug units; 
E6 is the restoring moment due to the elasticity of the wires and to the actual change of angular 
position of the model, and the term h ^ is the moment due to wind and frictional damping. 

It is the purpose of the following experiments to determine the value of 6 under various con- 
ditions. 

Before the value of 6 can be determined the differential equation of motion must first be 
solved. This equation is of the type of second order with constant coefficients, because, after 
reducing, 

(ID 

where _ 

bJ£, and C=f 

The solution of this type of equation is known to be of the form 

,-.(-!* V!R)< 

In order that the motion may be oscillatory the expression under the radical must be negative 
Then, 



The period of the oscillation is 



■ < 2 X« 



2ir 



and the time to damp the motion to^times the original amplitude is ^j^*?. 

T 

In these experiments the time to damp the amplitude to one half of the initial displace- 
ment was observed. Therefore 

, _ B-t l-g-t , 2 Log e 2Xl m 1.386 7 W 

If I m and t are known & can be calculated. This & contains three parts, the mechanical damping 
& 0 , the damping due to the wind on the apparatus arid thejlamping on the model itself, 6 m . 
To find the friction damping on the apparatus it is allowed to oscillate with no wind blowing, 
Next the wind damping on apparatus is found by allowing the apparatus to oscillate in the 
wind with the model removed, which gives + is obtained by subtraction. Finally the 
damping on the model itself is obtained by allowing the model to oscillate in the wind, which 
gives &=& 0 + &« + & m ; by subtracting (&.+&„) from &, & m is determined. 

The moment of inertia of the system is calculated from direct measurements of the modulus 
of torsion and the period of oscillation without any wind. This neglects the effect of friction 
on the period, which effect, however, is certainly less than 0.1%. If this factor be neglected 

„ 2tt 2v , I T*E 
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The direct measurement of E is performed by applying weights at a known distance from the 
center of rotation and observing the angular displacement. The value of E found in this case 
is .765 ft. lbs./radian. The moments of inertia of the entire oscillating system for different 
positions of model are calculated and tabulated in Table II. The -weight of the flat plate and 
the tailpiece are very nearly equal, so that the moments of inertia for same position of the two 
surfaces are the same. These values as tabulated in Table II are used in calculating the If,. 

II. 

Determination of moments of inertia of apparatus and model (at different positions). 
[Ncrr.— J-dlstance between center of rotation and center of gravity of model.] 





Time for 






I 


10 oscula- 


Period T. 


IjSlngs-ft.* 




tions (sec.). 




15 Inches 


52 


5.20 


.525 




61.4 


5.14 


.514 


8 Inches 


51.2 


5.12 


.510 


4 TnffTiWf 


61 


5.10 


.505 


0 Inches 


51 


5.10 


.505 


Without modeL 


51 


6.10 


.605 



THE DAMPING ON THE APPARATUS. 

The values of 6» in these experiments are determined for different positions of the spindle. 
The results are tabulated in Table m and plotted in figure 4. It is seen that the damping on 
the apparatus follows very closely the law of linear variation with speed, and varies as the 
square of the distance Z, measured from center of rotation to spindle. 

THE OSCILLATOR TESTS. 

On the flat plate two trial runs were made before the final experiments, the results of which 
are here recorded. For the tailpiece two readings were taken for every different condition of 
testing; namely, for different speed and positions of the models. The average values are used 
in the computation of M q , the values of which are tabulated in Tables IV to VI. Since 
Mq is directly proportional to speed, a simple basis of comparison is obtained by dividing the 

damping coefficient by the wind speed. The values of y 3 are plotted against distance from 

the center of rotation in figures 5, 6, and 7. 

m. 



Determination of damping on apparatus. 
[Now.— J — distance from center of rotation to spindle.] 
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U 


ft. 


30 miles per hour. 


f 


ft. + &. 


ft. 


15 Inches. 
12 Inches. 

i Inches. 

Oinch... 


375 
339 
283 
255 


.00193 
.0021 
.0025 
.00274 


ar 

230 
210 


.00256 
.00281 
.00305 
.00310 


.00093 
.00074 
.00055 
.00045 








20 miles per hoar. 


15 Inches. 
12 Inches. 

8 Inches. 

Oinch... 






278 
272 
243 
235 


.00256 
.00258 
.00258 
.00298 


.00063 
.00048 
.00038 
.00024 








10 miles per hour. 


16 inches. 
12 Inches. 

8 Inches. 

Oinch... 
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302 
258 
241 




.00222 
.00232 
.00270 
.00290 


.00029 
.00022 
.00020 
.00016 
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IV. 

Experimental damping coefficients for flat plate. 
[ I— 14.75 inches; from center of rotatlon.to center of pressure.] 
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6. 


, *• 


6, + 6. 
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30 


.0015 


.00093 


.0024 


19.8 


20 


.0015 


.00063 


.0021 


81. S 


10 


.0015 


.00029 


.0018 


58.9 
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5 - 4. + ft. 
+ 6m 


ft.-6-fto 
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ft. 
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30 
20 
10 


.0360 
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.0134 


.0336 
.0220 
.0116 


.00112 
.00110 
.00116 



{1-12 Inches.] 
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6.+0. 
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.0016 


.00071 
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20 


.0015 


.00048 




45.7 
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.0015 


.00022 
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ft-ft.+ft.+5. 


ft.-i-6.-ft. 


bm 
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30 
20 
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.02384 
.01568 
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.0071 


.00072 
.00068 
.00074 



[!- 8 inches.] 
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6. 
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.00470 
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.00030 



[(— 0; center ol rotation and center Qf gravity coincide.] 
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Fia. 1.— Lift on flat plate at vartoas speeds. 
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Fio. 2*— Lift of tall aur&ce, elevator at 0°. 



• 0 *8' 16' 
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Fio. S.— Lift ol tall surlaee, eleva- 
tor at 30°. 
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0 2 4 6 8 10 12 14 16 
I ■ Distance from center of mf of ion to spindle 

Fio. i.— Damping of apparatus for different positions of spindle and 
speeds. 
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Distance from CJ 3 to center of rotation in inches 

Fio. 6.— Damping on Oat plate. 
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O 2 4 6 8 10 12 14 16 
Distance from CP to center of rotation in inches 



Fia. 6.— Damping on tall surface, elevator at 0\ 



Flo. 7.— Damping on tall surface, elevator at 30*. 
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V. VI. 

Experimental damping coefficients of tailpiece. Experimental damping coefficients of tailpiece. 

[Elevator setting 0*. Wind speeds 20 and 30 miles per hour. Frietlona [Elevator setting— 30*. Wind speeds, 30 and 20 miles per hour. Frio- 
damping- 6,- .000778. Distance from C. R. to C. P.— 14.7 inches.] tional damping- .000927. Distance from C. R. to C. P.- IS inches.] 



V 


■ = : so 


20 


». 


.000930 


.000630 




.001708 


.001408 




85 


46 
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• .02065 


.01670 




.01894 


.01428 




.00063 


.00071 
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V 


30 


20 


b. 


.00093 


.00063 


&.+& 


.00186 


.00166 




29 


45 


b 


.0259 


.0161 


i- 


' .0240 


.0145 


y 


.00080 


.00078 



[ Elevator setting, 0". Wind speeds. 30 and 20 miles per hour. Frfctional Elevator setttng-30\ Wind speeds, 30 and 20 miles per hour. Frio- 
damping- &o- .000778. Distance from O. R. to C. P.— U.7 inches. J tional damping- .000927. Distance from C. R. to CP.— 12 Inches.) 
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[Elevator setting, 0*. Wind speeds, 30 and 20 miles per hour. Frtatlonal [Elevator settlng-30". Wind speeds, 30 and 20 miles per hour. Frie- 
damping- 60- .000778. Distance from C. R. to 0. P.— 7.7 Inches.! tional damping- .000927. Distance from C. R. to O.r.—S inches.) 
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lEIevator setting, 0% Wind speeds, 30 and 20 miles per hour. Frfctional ^Elevator setting— 30*. Wind speeds, 30 and 20 miles per hour. Frie- 
damping- .000778. Distance between O. R. and O. P.— 0 Inch.] tional damping- .000927. Distance from O. R.toO. P.— Olnch.) 
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PART ffl. 

COMPARISON OF CALCULATED AND EXPERIMENTAL DAMPING COEFFICIENTS. 

In order to facilitate the comparison of the calculated and experimental values they have 
been tabulated side by side. In place of giving the actual values of the damping coefficients 
in this tabulation, the average ratio of damping coefficient to speed has been used, the assumption 
being made that- the actual deviations from strict proportionality to wind speed are smaller 
than the experimental errors, and that the average is therefore nearer to the true ratio for all 
speeds than any particular observed value. 



FLAT PLATE. 
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TAILPIECE. 







Elevator at 0*. 




Elevator at -30*. 
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mental. 
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The check between the experimental and calculated values is quite remarkably good. On 
the average, the experimental values are a little smaller than those obtained by calculation, but 
the difference seldom exceeds the probable experimental error. 

There is a remarkable difference between these results and those obtained in another series 
of experiments made at Massachusetts Institute of Technology in 1917 1 , where the conclusion 
was that the damping due to the tail surfaces on a complete airplane model is 50 per cent more 
than the calculated amount. The difference can only be attributed to a systematic error in 
one set of experiments or to some fundamental difference between the conditions of the two 
tests. If any systematic error exists it is more likely to be in the first set of experiments than in 
those described in this report, as the new apparatus is a decided improvement over the oscillator 
originally employed at the Massachusetts Institute of Technology. It is improbable, however, 
that any of these experiments would permit of an error so large as the difference between the 
results of the two sets. As for differences in surrounding conditions, the only important one 
is the difference between a complete model and an isolated plate. The presence of the wings 

i AnlnvesUgaaon of the elements which contribute to statical and dynamical stability, tj A. EJfmtn, E. P. Warner, and G. SL DenBnger: Third 
k nniial Beport National Advisory Committee for Aeronautics. 
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undoubtedly affects the damping action of the tail surfaces, because of the reduction of aix 
speed behind the wing cell, the down wash, and the increased turbulence in the stream, and this 
effect may be important. Its magnitude can hot, however, be determined in general, even by 
an approximate rule, until much more extensive tests have been made. 

It is difficult even to predict from the data now available the nature of the effect of inter- 
ference on damping. The reduction in air speed would naturally be expected to reduce the 
effectiveness of the tail in respect of its contribution to dynamic, as well as to static, stability, 
but the question of downwash is more difficult. In so far as downwash is a function of angle of 
attack alone, it should not affect damping if the axis of oscillation of the model passes close to 
the center of pressure of the wing, so that the true angle of the wing to the air is not affected by 
rotation'. It should be remembered, however, that various points along the chord of a wing so 
mounted are moving in different directions, and the result is much the same as that of changing 
the camber of the wing section, the effective curvature being deeper when the pitching rota- 
tion is positive, less when it is negative. The lift coefficient and downwash at a given angle 
therefore appear to depend to some extent on the angular velocity, and their variation with 
angular velocity must produce a secondary effect, probably very small, on the damping due to 
the tail. It has recently been pointed out by Cowley and Levy that the time lag before the 
downwash reaches the tail may have an important effect on the damping coefficient. 

In short, while these experiments leave much to be settled by future research, they do at 
least justify the methods hitherto employed for an approximate calculation of damping, and 
they show that the damping directly due to rotation is negligible by comparison with that due 
to translational motion of surfaces as far away from the axis of oscillation as are the tail sur- 
faces on an airplane of conventional design. 



